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PHASE DETERMINATION OF A RADIATION WAVE FIELD 
Ftetdofth* Invention 

This N mwa;aikm wa 

Fix Iyengar- also a? pha 

about a mdmtioa wave field can ba used. As used In this specification the term 

^ A - - > ^ a , - \ m,^ ^ ^ N ^ 

in a wave like manner including but not limited to examples ouch aa X-rays,, 
visible Add and electrons. 

\ or i m ; o for thB measurement of the phase f 

inany aualiaatloaa In fandarneafal physios and as a basis for a number of 
- N \ \* of 

applications of phase measurement techniques include the fields of x-ray 
imaging, elecfaea rnicresoopy : optical ^ as wall as optical > w 

and x-ray phase tooeyaaehy 

Phase is typically measured using Interferometers of venous types. The key 
oafui-a a t mvwet ear tv to qualitatively m-osme the hh-*se of a 
w - e A w-. - » \*sed teehnkfuea retain significant importance 

If has been recocmsed that ncw-ime;Aaona.Anc techniques .may ne used to 
provide phase information A combe? of mam miasm mm etna approaches involve 
attempting fa solve a transport of inteaslty equation for a oca scan wave field. 
This equation relates the irrsdiaoce and abase of a paraxial monochromatic 
wave fa its longitudinal irradlance derivative and is described in PAR. league, 
"DetermiiasiiO Phase Peinevsi, A Green's function Solution" J, Opt Soc. Am. 
73 '1434-1441 (1984). The abide "Phase imaging by the transport of intensify 
equation" by hi. Sbeiei. Opt.. Comm. 4S 6-10 (1984), describes an approach 
eased an the transport of intensity equation by which phase struoiuee can be 
rendered visible by the use of defocua and digital subtraction of mieositv cats 



obtained at various defaousad distances. This approach only provides for phase 
^ ,v > v ^ \ - \ a s ^ 

approach based on sowing the transport of intensity equation is disclosed in I.E.. 
Gureyev, OA. Nugent, D. Paganin and A. Roberts, "Rapid phase retrieval none 
$ a Fast F< s i O , h RavOmbp on lew o 

daedca Te , s n K A 

Nugent, "Rapid quamriafoe phase rmagsng uasng the transport of intensity 
3 i 133 m 34 

a Horn held to Po recovered from two sioseiy spaced intensity measurements 
to - bsam has an arMfrarv v - - 

dlsihbufion limited by rectangular aperture. Whilst this approach can be used for 

non-uniform intensity distributions the extent of the non uniformity is Notice and 
s ^ s v u&onal complexity, Como , * 

not able to cope with non unites; uses introduced by some sample absorption 
15 prof m > i some 

aisaaniyaap :abte t snt wave flsHfe. 

The article IRA Nugent, T.E -Gureyev, >D.f\ Ccwksen. D. Pagans and Z. Barnea 
s . ^ m using hard X«ray$ H (1996) 77 Rhys. Rev. Lett, 290 1- 

20 2984 Is also based on a solution to the fmnsport of Intensity equation. Again the 
technique described cannot be applied to a non-iifWornt intensity distribution, 

m ^ \ * m a solution to the transport of Intensity equation 
irmbed to a requirement of uniformity are described in T f. Gumyev, K.A.. 
2S Nugent, A. Roberts These retrieval with the franspart-of- intensity equation; 
matrix solution vain the use of ZernNe polynomials^ J. Opt. Soc, Am, A Vol 12, 
1932-1941 (1999) and I.E. Gareyev, A, Roberts and K A. Nugent "Partially 
coherent fields, the transport- (A intensify equation, and phase uniqueness", J. 
Opt Rom Am. A Vol 12, 1:942-1946 (1 995). 



described m I.E. Guryev and K A. Nugent "Phase retrieval with the transport- 
ofdmensity equation. II. Orthogonal aeries solution for nonuniform m ^ 
J, Opt Soc Am. A Vol 13 NN hi )u f nn war , s ^ 

5 s vt\ in 
implementation. In many applications - compiexSy makes the technique 
Impractical, 

^.■-U'sviCtsiuo m i.5 5 e ii!v»!!wO!> 

the measurement at phase, in combination with a caeci measurement at 
intensity a measurement .of phase allows the phase ami imvnaiw ai -any othur 
plane In the radiation wave said to ae determines using known techniques. The 
Invention - provides the basis for a lumber of moo- , , s m:m> 

15 

in accordance with a first aspect of the invention them la provided e method of 
of the phase of a radiation wave field tuning the 

•stops of 

■(9) p <H ^ » measure of the - vm m ^ . m a . nrw u . 

20 of said radiation wave field aver a selected surf&oe extending genemlly across 
tiio wave a em 

fd) producing a representative measure of Intensity of said radiation 
wave field eve? sold selected surface; 
fa) v °* - s s 

m erst integral transform representation and applying to said first integral transform 
representation a first filter corresponding to the inversion of a first, differential 

I change asm am to produce a 
first modided Integral transform rapmsentatlon: 

(d) applying an inverse of said first integral transform to sard first 

:.w modlhed integral transform representation to produce an urmansfmmed 
ronresO'rOai.;or;; 



4 

go) applying a correction based on said measure ef Intensity over said 

selected surface to sa;d untransformed mgresentatlon; 

(i) transforming the eerrecfed antransformed representation ta produce a 

second ^ v transfara^ representation and applying to said second rnfoara! 
transform representation u seeena filter corresponding P> the inversion of a 
seaaad - i fleeted in Ida corrected untransfon-ned 

educe a second ^ m ; v ^ s < > . - N 
(g) applying an inverse of said second nter ras transform d: a ml s< cos d 
wo'Mzd integral transform representation to ax »e o 

said radiation wave field across said selected plane. 

; accoidr oe with a second aspect of the invention teem -a ivovided an 
apparatus for c|uanf»a detonation of the phage of a radiation wave field 
including 

(a) means to produce a repmsentath/e measure of the rale of change of 
m Lie o.av i del o\ > * seler m soda 5 

extending generally across the direction of propagation. 
{b) means to produce a representative measure ef intensity of sad 

radiation wave field over said selected surface; 
(o) pmcesslnr aa s to sequentially 

(I) transform said measure of rate of change of intensity to 

(II) apply fa sain first integral transferal representation a first filter 

c i the wersk v of a 1 si 

said as da a { aensip to 

produce a first modified integral transform representation; 
did apply an inverse or said - integral namdsm to said first 
oedlfied $ sfon epres* Yfatio < > or c« an 

unuanrdorrned represenmtiora 
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(IV) apply a correction eased on said measure of intensity over said 
selected surface la said umransformed representation; 

(V) transform the corrected uafransformed representation to 
produce a second integral transform representation; 

(VI) apply to said second integral transform representation a 
m nd ftltet correspond 

differential operator reflected In ihe corrected ardransforrned 
mpmssxmm m patter* second modified integral transform 
representation; and 
N\\) appiy an inverse, of said second integral transform to sad 
second modW misgml transform r-epraaentatsoa to m xmce a 
maasore of phase of said radiation -wave- field mom m® 
detected plane. 

The selected surface can take any form thai extends across the direction of 
propagation of trie radiaHori Including planar, psdavpheneat and padmyllmrfrlcal 
surfaces:. 

The first and second integral transforms can be of any suitable type ami delude 
approximations employed for oomputattonal convenience, speed or efficiency. 

The first and second integral transforms are preferably produced using a Fourier 
transform. More preferably, the transform is a Fast Fourier transform. The 
methec and apparatus of this mvervsc^ provide for determination of prase of a 
radiation wave field In a manner that Is computationally significantly less 
complex than prior an approaches. This results if. sroodcantiy lower 
cornpidasan times in some examples computation times improved by many 
v ^ s - -m * ' mu-m- mac- 



The first and second difieremial operators are preferably second order 
differ* * > s natation of the method and 
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apparatus the first titer la substantially the same as the second filter. It Is funtaer 
preferred that at teas! one of the first and second niters Includes a correction for 
^ ^ el Intensity 

o e f< m of tt« went >n the first filler < > N < « 

s « st higher frequencies of the first Integral Ira pm 

of the Invention the second filter can comprise seh\. . suppressing second 
higher frequencies of said second integral transom repejsemanon, 

•The correction based on me measure of intensity over a selected surface can be 
a nil correction where the Intensity mn^m are a a predetermined 

selected amount 



ne measure of the rata of change 

< i d surface are produced from sneasi. env^c, w m-- 
Intensity distribution over at feast two surfeces emending across the wave field 
and spaced apart along the direction of propagation of the radiation, in another 
form of the invention the mpmsentative measure of rate of change intensity in 
the direction of radiation propagation is produced by obtaining a first 
representative measurement oyer a measurement < \, <te« ga s? \ 
>n fo radiation s N v 

representative measurement over sakl tmmurement surface for radiation of a 
second different energy, in the case of X-ray radiation, for example, the change 
In radiation energy can be achieved by changing the X-ray target or by suitable 
filtering. 

The selected surface for which measurements of Intensity and rate of change of 
Intensity are produced Is preferably located between hvo of the spaced apart 
surfaces over which intensity distribution Is measured. 



P€TaUI9MKK*4<? 



In bm prfemred s at the snve-t^on v Elected sedans- - spaced apart 
surfaces are planar. . it is further preferred that the planes are generally 
perpendicular to the average direction of propagation of inn radiation 

s The raeihco and apparatus of this invention can he at least partly implemented 
us-np a suitably programmed computer, in canicular the processing means 
preferably comprises a suitably programmed computer and the steps of the 
method are preferably performed ue ^ . . ^ ^ . * In 
such forms cf the invention Intensity input information may take the form of 

•o deatfeed images or data containing Information from such Images, in other 
implementations of the invention a dedicated Fast Fourier transform chip can he 
employed as at least part of the processing means. 

TK a ^ move iseasuraof rate of c&anpe o \e - 

a; by sefetraoiion of representative measurements respectively made at ioostlons 
over one spaced apart surfaces. In the preferred form of the invention the 
representative measures of intensity and rate of change of intensity are obtained 
by sampling measurements at selected locations over the surface. Preferably 
inesampiinc and. measurements me tm$& #'J©es8©ns defining a regular array 

2.0 over the surface. This can be readily achieved for example by using a CC D 
> n ^ M^deteoto 

in the preferred method the direction of propagation of the radiation wave field Is 
selected to he the • direction at a Cartesian co-ordinate system and x and y 
:>5 components of phase are produced separately. 

his Ca s > or J mate syst m do i is the direction of 

propagation of the radiation, the preferred filters are of the form 



<Af+A i 
W a;.,y:)m,^y 



o _ 

e 

5 k*, k y am the Fourier variables conjugate to x y and 

a Iwammimtm^mnmdb^ mim mthK \ x - 

Is equal to zero for a no noise case. 

o ate f o mgc the 
?<5 i * of the average wave number ot 
< ^o into the Foiuk » a — 

Tne representative measure of Intensity ov»rthe spaced apab surfaces can be 
obtained by imaging of that surface through m appropriate imaging syaiem 
l5 >st 5 s nansi^ formats i ma be k\ < * than 

measured at the surface. 

The method of this invention thus provides for the quantitative and decoupled 
determination of phase and intensity of a radiation wave field at any surface 
across the direction of propagation of the radiate. From this phase and 
intensity deternsrrason it is possibie to calculate the phase and intensity at any 
other surface along the direction of propagation. Accordingly, the invention 
provides the basis hv a number of measurement techniques. 

25 in a fuUnet aspect of the invention there is provrdeP a rotthod of mugtng an 
obiect including tne steps of 

(a) exposing the object to a radiation wave field from a source, 
<b) producing a representative measure of the rate of change of 
intensity over a selected surface extending generally across the 



wave field on the < > 

radiation; 

(c) producing a representative measure of intensity of said radiation 
wave Held ove; said selected surface: 

g - aid measure of s c to 

s leg-al eansfo eproseno sppiysng 1 

said first Integral transform representation a first filter 
(Responding to -the Inversion of a - diff&remiai operator 
reflected m mM vtmmm of rate of change of intensity to 
p«>duGfc-« fffst r?io<»^.-S^gi#i:tran$fo? n repu sanitation 

(e) applying an inverse of said first integral transform to said first 
modified integral transform representation to produce an 
uniransformed representation; 

(f) applying a corrects on based on said measure of intensity over 
ssk! seteeted surface to said untransforrred reps jsantafio - 

v \ - - g the corrected untra isi a med \ - v to 
produce a. second integral tra 

apely ; ng to sold second Integral transform representation a 
second filter corresponding to the inversion of a second 
differential operator reflected in the corrected urararcrformed 
representation to produce a second modified integral transform 
representation; 

(0) applying ao Inverse of said second integral transform to said 

odnoo a 

measure of phase of aa?d rad«atlon wave Held across sad 
selected plane. 

In a still further aspect of the Invention there Is provided an apparatus for 
maejlng an adject hi dreg 

(a) a source to mediate the object with a radiation wave field, 
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ftp means to produce a representative measure of the rate of 
chance of intensity of said radiation wave neid ever a selected 
surface extending generally across the oeve held; 

tc; means to produce a representative measme of Intensity ot cam 
radiation wave nod over said seieoted surface: 

gag processing means to s 



(0 



of rale of change of Intensify to 
produce* firstihfegm ire 
(H) apply to said first 'mfegmf transform reprosentatd-n a 

mm mmzpommg to inversion of a first differentia! 
opemtor refteefed in said measure of rate of thmm of 
intensity to produce a - rrmoifiee Integral Pansi'om-i 
representation; 

uppW an inverse ofsakl first integral transform to said 
modified integral transform representation to produce an 



(Hi) 



dV) apply a odrre#fart based on sad measure of intensity 
over said selected surface to smd unmaoo\nsnvd 



(V) transform the corrected urttmnsfofroed representation to 

produoe a second ire 
(VP apply i.o said second Integral transform representation a 
second Met corresponding to the r^ersaon a second 
operam ! > cameUed 
reservation to praduce a second 
si transform representation: and 
rviP appiy an inverse of said second integrai transform to safci 
second modified integrai transform representation 
produce a measure of phase of said radiation wave neid 
across said seiecmd plane. 



to 



The radiation used to Irradiate the object can be a planar wave bald or »pherical 
wave bald or an arbitrary wave field, if It Is desired to reproduce the chase in 
the obiect plane the phase wave « determined by tbe abova meihed and 
s r ibis 5 >ackr. opagaa. -d inc eld us la u 

Tbe method and apparatus of Imaging substantially incorporates the 
ceiem-macm, of phase as disclosed m Nation la the first and second aspects of 
the invention. The preferred aspects of the invention described m relation to 
nose aveects n we a a also applicable to tie method and apparatus of 
imaging. 

It Is ne^bie « some appMcafiens to use a -m* object to image plana distance 
eonespondlng to contact-imaging with asm propagation disbmce 

If M IN object an i\ eoonstiucted in \ v \ %m by back 
oroaagatiag the intensify and quantitative phase informal to numerically 
'reconstruct an image of the actual object phase and intensity Mructu* 

th other terms of the method more than two image plane intensity distribute 
ts can be made to obtain a belter estimate et , - ^ * 
tensity w intern va i this case one o both of s< >fc « r 

< u , m ^distances is changed an. a ^ ^ — 

> natk * or. s recanted , 

measurement Is made Tan measurements provide eats to which e function 
^ m av^c^mga - He-wy. 

The method of imaging an object has particular application to point projection 
microscopy using X-rays visible light or electrons. 

■ In another aspect this invent provides a method of phase amplitude imaging 
including the steps of 



(a) irradiating an object with a radiation wave field; 

(b) focussing radiation from the object through an imaging systen 



(0) - !i ^ve measure as smenssly 
distribution of radiation over said Imaging surface at a ftrst 
focus of the m&gmg system; 

(d) introducing a change in focus of ton image on sak! imaging 
surface through the imaging system; 

(e) )ro< c \g a second sprearm 
distribution over said : mu ng mmmce; and 

{!) using said first and second mpm* ^ -->\ to 

t^m^^m measure of rate of change of intensity over a 
selected surface extending across toe wave field; 
transforming said measure of rata of change of inreroa^ to 
produce a first iofegraftransfonm representation ami applying 
to sast first integral fransfom> ^ . a ^ a brat niter 
corresponding to the inversion of a first differentia; operator 
reflected in said measure of rate of roange of infeossfy to 
orsame a first modified Integra! transform representation; 
(h) anplying an .averse of said s k ; rsmpmi Worm ta , ^ 

anfransforraao representation; 

(1) applying a correction based on said measure of intensity over 
said selaetad surface to said untransformed representation; 

(g transforming the corrected unfransforrned representation to 
produce a second integral transform representation and 

jppjyrng to said seer gra! transfc < ... tat a 

second alter corresponding to the inversion of a second 
differential operator reflected in me somsmno >^m> .uonaod 



S3 

representation to produce a second modified integral transform 
representation; 

(K) .applying an Inverse of said second integral transform to said 
second modified integral transform representation to produce a 
measure of phase s < said mdlat )r warn field across sard 
Reefed piarm 

> ■* a psc of m Invention here is provided an ? p? at s fo 
phase amplitude imaging of W object including 

a radiation mm field sour< e k ao r sa d on 
an imaging system to focus rae % 
imaging surf&oe ex 
object; 

means fo p oduco a ep esentatne o 
, ^ x - s^d imaging surface; 

said imaging system Induding selectively operable means to 
adjust said focus of said wfafon to said imagine surface to at leant a 
first focus and a second focus; 
processing means to: 

* v\ . s\ n ^e x - a 

representative measure of rate of change of Intensity over a 
as morn too vmne - from 
mptesentaffve measures of radiation intensity over nam Image 
surface at sale first toons and sale second focus; 

(II) transform said measure of rate of change of intensity to produce 
s first integral transform representation; 

(III) apply to said test integral transform representation a first filter 
corresponding to the inversion of a first differential operator 
reflected in said measure of rate of change of intensify to 
produce a first modified internal tmmdono representation; 



(iy) apply an inverse of said first Integral Imnsfona lo ,m : < m ; g 
modified integral transform representation to produce an 
untrancmrmeo representation; 

(v) apply a correction based on said Pleasure of intensity over said 
Selected u 

(vn , ae^eP-v* - to 

produce a second integral transfb? 

{vip apply P> sale second integral transform representation a second 
filter corresponding to flie inversion of a second differential 
operator refected in tb@ cop<, , ^ ^ , 
to pmdupe a second modified integral transform representation; 
and 

(vi) apply an inverse of said e^ s \ c . a mo m . 
second moeifed iniepal traosform repfe^a^ to v s 
mmmm of phase of said radiation wave field across said 
selected plane, 



Preferably, the numerical aperture of the sradiahng wave held Is smaller than 
the numencai aperture of the imaging systfem. 

Preferably, the Imaging surface Is a detector; The defector h of any suitable 
ferns sum as for esampis a CCD camera. 



)af;ai resolution Is minimised. In scree applications more loan Pan -ma 
s obtained to obtain a Defter estimate of the rate of change of intensity. 



fee N <* fiw pnasse ^ , \ 

aaofparates the determination of phase as disclosed sn relation m the ttrsr ana 
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v * " ^ ^ < N <N >e0 aspects O the *mo\m x 

* aspect; ~- 1 m- x moth xl 

and anaemias at imaging. 

in a preferred application the method is used for quantitative phase amplitude 
microscopy, in this case the r , i 1 

m tea preferred form of the ifmt&m the sudaco >s me?ambiy planar. 

The invention will now he further described by ws\ o! example >ni\\ with 
mecreme \ \ J* > , v» m wbM%: 

?iQ%sm 1 js a schematic illustration of an arrange.rm a for i 
1 r aa- \ N an object is illuminated with (a* piano 



Figema 3 fa} to fit are simulated Images siiuiUrMsng phase aamsaioahon ror 
) ic e *v 

m - - a ^mo{*ss| are a safes of images 
hack propagation to another image plane; 

Figaro S is a schematic representation of an arrangement tor polar projection 
necroscopy using the method of has invention; 

Figure 6 is a schematic illustration of an arrangement for quantitative phase 
25 amplitude microscopy using the method of this Invention; 

Figure ? is a sohamatia emvann of an exemplary system for quanfrlatlve phase 
amplitude microscopy according to this invention; 

Figure aha to hid show Intensity images and phase images obtained using the 
system shown ;n Figure ?; 
50 Fiqur** a ;s a graph showing a comparison of measured a\ expected en a; 
profiles for the fibre at Example 3 



wo emtez?. 
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Figure 111 is a schematic drawing cf an exemplary system for three dimensions! 
optical pha w accr pa I weeflon 

Figure 11 is a schematic enlargement of part of Figure 10; 

styles tomograp o slice throng dme \ t m or in 
s Example 4; and 

* ' o" " ^ - m o s , - as - - 

known refraetore index distribution according i:o Example 4. 

Figures 1{a) and mi show a schematic aaangen-erh far phase deteoninahon in 
io vac* hv - with this lnverto : - where an object is illuminated by pianenvave 
radiation 2 or point sounm. rad&Son 2 to produce t f « px n-jeas s 3 

V each point in space, an op Pen ace a 

oaase. Intensity Is a measure of i&e amount of energy th \ \; p, , *«, v 
hi o^ . o - , measure of the direction of the energy flow. 

internet-/ may ha measured directly, lor example by recording an image on flier 
I * \, v s s isvasnred using interference vetn a "reference beams In 
_ a \ ^ \ method gives a ncn^merferometn'c v , i 

intensity can he measured over hvo parallel planes A, B extending across the 
direction ol propagation oi the wave Paid on the sloe remote trooi the Incident 
radiation > 

25 

'The preswd invention determines phase by providing a solution to the transport- 
oplrhenelfy agnation: 



(1) 



V «0Ptf } 



n 

where I is the intensity in the plane, the gradient operate m tne plane is denoted 
v . k is the wave emceer of the radiation, and ditdz Is the intensity derivative 
or rate of change of intensity. Mete that ®ldz is estimated from five difference of 
the measurements in the pianos A & B shown In Figure 1, white the intensity I is 
s given b\ dx sve age at I 

In ardor to obtain a solution to equation 1 the « 

m i ns \^ v . N * necorn^: 

(3) Y, •(V x A)~~mj\ 
mh ng i m a* a standard de MkyV^s ?/jhisn a> i e * ae i 

15 

(41 9?d--4»Jl 

« ^ x vai-ohv -\ sjonal Upiacten acting over the surface of the 

:mage. This equation has tne foO:av;na symbolic solution:. 

20 

f the gradloal operate s apt ) sides c uaho becomes 



! ne do >g equation N a tie fund - a 3 ;V v<s ' s t a-h mod Into 
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Dividing both sides by I then yields; 

Taking the iwo dimensions; divergence V * of both ssdes of (8) : and again 
making use of the standard identity \ *\ - v then (8) becomes; 

This equation has the following symbol mMon: 

(10) 4 - ~*v.{^ ♦[/•• $ V:V,; J a,/J. 

in order to implement a practical solution to equation (10k the foliovdng formulae 
are required A suitably well-behaved function f(x,y) may be written In tive form 
^ - - togral: 

(11) my) - J J /<*, Ay^dk^ 

The function }\k^ is called the Tourtof transform* of fix.yj. 
Thexdenvauve ef (11) yields; 

(12) | ffe.y) I f | n x h \ ^ I 



a 

Henc« < >r transform of 

fboy) \ o \ by * < Stated differently, - fF %? ( where F denotes Fourier 

OX 

f denotes 
considerations apply to ~j~f{x,$ 

5 if the Laptacian vUjL, |™ of {11} is obtained and similar reasoning 
applied it follows that V* t where kf - *; 4 *; . Thus: 



(13) V r » »F*U% 'F, 4 85 **"V< T" 855 *F%F - 
Sj? ay 



jo Hera, F denotes Founer trah9&nnaliot t , F' 1 denotes 
1 5 Equations ( 1 3) can be used to rewrite equation (.10) In the form 



(14) #-#«+^» 



In praobce division by Intensity is only performed If that intensity Is greater than 
a certain fi r N 5 rvalue) 



Instead multiplication ay sera takes Place at this point. This amounts to taking 
the Gar >alvak Pth .... 

kn to 
Incorporate some physical constants into the phase recovery algokihm given In 
Equation (14) reiaknu to the experimental setup In use to quantify the variables 
k X; k s This can ba done by rewriting equation (14) in the following form suitable 
loi *h*ntv*v ^ ^ -9 a cast Fourier transform. 



10 #^4-^ 



h v ™ ™ J index the- frequent ampomtM of F{1, - /..) 
intensity derivaiive , s >. obtained by subtracting two kceqes h and L 

N VV P 

is uang ma i a 1 that the Fourier spas® step size is giver? ay 



v/her-v the Image is sn si x At array of pixels of shoe :\v. Thus in addition to 
measnrlng the three intensity distributions k Is necessary to know the pixel size 
a.? < defaces distance a> and wavelength A in order ta reako a quantitative 
phase measurement All of these quantities can be readily determined: the pixel 
sen he determined daeotk - CC to vmetr> (in 

the t em imagsng by a.- i the 

< \%\ a sP sc,aie< mse of an aging sysien iefoeus 

PsP \s N > v s N - s v c%\ and e - < 5 3 > o1 the 
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llkawnahon can Da determined ostrser ay monoehrornating the bicideat field or by 
analysing the spectra; dicfncuhen of tne mdoumn using existing spectroscopic 
PiethOdS. 

5 An example cf the phase-retrieval method implementing the solution of equation 
(14) can be representee by the flowchart shown in Figure 2. As shown in Figure 
\ <ni of phase o! a 

a set o muoomv measurements hp I oxer the two spaced apart planes A and B. 

A measurement of central intensity kx./i In a selected plane parallel to end 
m midway between the planes A and B Is also obtained. The Intensity 

^ v ittribtmtd over a cieHnedl array on euA: <,t Ao two planes A 

and 8 and the respective value* subtracted to produce a measure of the 
dive Tea; vak $ ultkdk vvaw 

m. - I it* ~m x- k mto two eompcm Awe \1 
w A performed to produce the respe x s x 

domain. A filter Is then applied to the Fourier don\ to 
don of a d fen p c the 

untransferrned representation. The eifterentAI operator Is represented by o;AA 

20 traeslaeo > Ten performed on each of the X ami \ ponems 

spatial domain value from which the differential operator has been removed. A 
dMsfcm by the central intensity kx.y) obtained by averaging the Intensity 
measurements ever planes A and 8 Is then performed r? the Intensify level is 
t , ^ \ v » me &i rani data is aga 

2.5 IrensfoericP ana nmleelleu by the same id no to aga a correspond to the 
> r. < * - , A s < ) he: s The 

> ^ , v ^ s ^ v \ - and 

A'vTx for the y component . The resultant components are again Inverse 

30 
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It c Ihi according th« bon can 

Koeeed from any suitable tope envative o 

extending v the 

propagation direction and the Intensity over that same surface- As will bo 

Usined h\ a vsnat :>f ways 
and the method implemented to yield phase of the radiation wave field. 

R events no conation nod with: 



oho v) denotes the recovered phase, 

F denotes Fourier transformation, and F' 1 denotes inverse Fourier 
hanse. s \ v 

l(x v y) le the Intensity distribution over the piane ef Interest, 



p ess 




20 where: 



(x,y) am Cartesian coordinates over the plane of interest 
(wmy) are the Former vansoies conjugate to (xw) 
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k « 2e/ I is the average e/avenumbsr of the radiation, 

, ; l ie toe aveeage wavelength of toe radlailom 
abas is the estimate for the longitudinal intensity derivative, 
e is the regatauzaiioa parameter used to stabilize the algorithm when 
s no so is en 

As given above, the solution to the transport of inters a equation (1) asm mea j 
, ; m an is, there am m > lh< Mot 

system used to obtain the intensity data which b ' \ - 0 Of 
to ^ uvaeet IN npesbx 

by a set of numbers: 

(16) 4,4,4^., 

13 watch are termed aberration coefficients, 

if intensity data ware taken on an aapebeai: tnatmrnent waose aaaerfeohoos 

;a«n set of mmm abem 
ii would be deseahfe if be filters fysgma: and bt ^ ,m present in (IS) 
20 eoufd be mpfaeed by modified filters which expilastly depend upon the aberration 
coefficients: 

(17) v,X«o\A. c»,.t,..fj,^,...> and ;dv^ v m -m m m » 

25 This would aiiow the imperfections of the Imaging system to be explicitly taken 
mats •> ^ ^ s sing aberrated 

imaging systems.. Far the aaealal ease of a non-absorbing phase object la a 
radiation wave field of uniform Intensity with weak (I.e. much less than 
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2* radians) phase sanations the appropriate modified lifers load to the following 

da) fay) - r _ t 

3 wham: 

tetmdM i» the aberrated intensity measured at detenus 

A, 

t •> <* 0 PN 

imaging system. 

to 

If a IPter Is defined: 

(19) Q^i^ t 4>^4,-.) s 

J 

is 

Then (18) becomes: 

(20) ,^.v^F^-F-V^U_. v » m 

The lerm - N * change o k 

2d o msliy is s ema^ ^ same 

p a: a i as (1 C smiueaiis the s be dealt 

with by changing the f liter in the eenenP method described above. The x and y 
component Piters .Pp and CP, am aivea by 
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A simulation was conducted in accordance with the arrangement shown In 
5 \ , N ^ N ^ N v x x s s * * ^ N the 

operation of the method on simulated nolsee'rea data. Diffraction paiiarnc are 
calculated using the mrmuiarmpectryrm formalism, an orthodox procedure. 
Figures 3(a) to 3(f) show images produced in the simulation. 

m \ c N 1.00 cm squam ami p - , ^ - < : ^ i a, 

x \ ^ * plane extending perp^icul-srh -con:; the n opagatlon 

direction of the radiation The wavelength at the light was taken to he 632.8nm, 
Tha iotencity e : a,e plena ? ~ g. which varies ham 0 to 1 ;n arbitrary on.m, ie 
shown in Figure 3(a). Within the area of nonzero illumination, the minimum 

- . v m Image cortex 
varies from 0 to a radians « shown in Bgure S(fe). 

$ to ptetm. negattvesiy and po< ^ ( $ hum 

28 t N >s m Figures 3{cj \ h,w* 

vt« of 180 and i?§ ah , 
induced phase canadas! is clearly visible in each of these images. The two 

c meaed images a fed to form tonsils givernvhioh a show 
^ Fig, ;r e 3(e). 

25 

Yosgas shown , . N sares of 

a asit> aiid ate of change ci up -sh across the ,-"0 plane were then 

oessod secern :n a com e m >d shown in 

5 3 f) Note 

30 hetFkgr I m c same greys eating that 
the recovered phase la quantitatively correct. 
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Figures 4 (a) to (h) shows a sanes of simulated images iimxaaAm phase 

, OA \ ^ ,*V Vv s All 

Images am 2o6 pixels x 256 pixels 1 cm x Aam >r dimensions, with the 
5 s a 532.3 en » a >asx the 

•> m a gneu ; - - N ; N \, ^ s 

:SDecta*e\ ^ s ^ at 

anoes m 20C of 
Irdenoation ham Figures 4 (as and (a; A Ae imeraay rnessumaieaA of Pkoums 
10 4 to), AD and tea Using the images of Figures 4 (c). id) ana (a) amy. tbt ph sse 
retrieval algorithm obtained the phase map given Irs Figure 4(f) for the phase of 

i a suance 200 mni Image* of Fig*, 
used to numerically back-propagate the field back fa - 5 gav« 

I a; r * 4 g - d (h) lot rbe back propagated intensity and \ 
;5 These are la excellent agreement w |th F Auras 4 (a) aria lb), thus demonsiraaaa 
me use ilau uu a x techni<|oasfartheqx s ^a s \ 
m>pillode ana abase of a tk se ovc which 

s e aAr Ibat I v I not 
restricted to free space; back-propagation can else be affected through a known 
I system 

As s x o\ s gua n a < m from a 

25 v \ s N N ^ " x ' ^ v\ 11, 
located at a distance d ;sc bam the scarce. The radiation passes though the 
object 1 A and Is allowed to propagate a AAA.u distance a- to ana of Ibe Image 
planes Ah b m wAcn tna xaemmy at Ac mash An a detected, in A detection 
A r- g a standard device such < x a or 

30 \ v a * v One 

or both of the distances a,,, and/or 4* is then changed so as to introduce 
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defotus Into the images and the intense distribution is measured ones again 
The case of corresponding to comsoeemtmng ,vpp -em propagation 

distance is Included as one possible measurement. 

5 The intensity data is then processed usmg 

recover the decoupled intensity and oha;>e information in tea sms^ng plane. 
Parameters, such as wavelength, pixel slam and defooas distances are inserted 
> age mas explained above, to yield qutm&ftv-: ^ovaat^n admit the 
magnitude of Urn chase shift In the Image plane 
;o I cr-nsir; oases a reconstruction of the obleesHnthe^ - * , * 

ihe downstream distraction plants U. I is des red In this . no r e m smm u t 
x N v - en obtained above can be used to ba i } \ 

\ — * the )b \ | ana thereby numerical ^ s 

she oduid object phase and intensity structure. This can be done using 

" to take more than two images m or des v>. omain a 
hot m m imonslPy derivative dl/dE, in which ease cn« or coth ot the 

distances d : ,, and/or d^ is altered enee agasn ax- 
so - 0s \ .^\\ the number of desired Images is acquired. A 
function can then be fined to this data from which dl/dx can be computet! and 
used sn the phase recover/ itepoOnm in nines os the oirstee subtraction ot two 
smages normally used. 

•25 Exjpgole . :j e Pmyppteyp. p ho so pp. v ires is . ; r ops eonpv 

FiguiB 6 schematically shows an arrangement for quantitative chase amplitude 
necroscopy, A sample Is illuminated using a source cf white light KShfer 
YPrmik or onlyfo ^ N * x scopes e light 
30 throt^h an object 16 and collected by a microscope Imaging system IT end 
relayed to a CCD camera IS or other digital imaging device having a planar 



2S 

imaging surface. Three Images are nnllsoled: en in-focus image, k and two 

v^^wi.v^ *a ^\ r-.Cc^ m - , 
such as $ tii se system IS to ed est the mlcrossr pe fc ms km b m di focus 
iMrpdiioed is usually quite small sc that degradation In spaUai msoiudaa is 
s Td raised, although vh c N mi ss by 

sample properties ana snspmg geometry sacn as magmtsmPom numerical 
aperfoa 



yVhca taCng the images the numerical aperture of id- condenser is chosen to 
10 belssstda ira numeiica) apertum of the objective doing used, if this is not the 
ease then cmfom image degradation WW occur* a8ho«gh foe precise namem 
which the condenser and objective ftunnerfcal apertures sPouih d^fer involves a 

/and spatial msoiution. with t 
lepeadincr N andfhe opticVused, 

15 

c v cted images I* arse ^m , . N . s N 
representative measure of rate of change of Intensity' amfomnty derivative). 
I \ \\ -mn e t> ^\ > ssry data of oolleeied imago \ 9 toe « ietd -c d- so h fo 
fx » - used to produce quantitative Mormstksn about u -e magnitude of 
20 the chase shift in ^ in-ege plane 

As la Example 2 for point projection, there may be oases in which it is desirable 
ia lake mere ma dm nag* * order to obtain a bettor estimate o ; itenslh 
derivative ^ ^ s s dec which 4Mz can 

subtraction of two scapes narmaliy used. 

It >ss i. - N \ to fota 5 

surface topography. The principle of operation is the same, however die optics 
so ve to b c k. o ^ e s to f o c v n s 

the process m identical. 
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For certain applications I can also be desirable to filter the light to a earttcelar 

wavelength) aifhoug ithfeis 

it warns equally welt with white light, 

5 

An expenawnial wpfienwctaiton is shown in Figure I' An Oknpaav S#-S0 

optica; microscope 20 was onutw 
. x - - andenser to pax , 

be able ib compare the results with extsilnr 3 j xi xx.aaan 
is and a set of cow-sip corrected Uplands objective v*e \ - — 

microscope, enabling mages to ba taken of the same field of view using bath 

phase retrieval and Nomarski DiC for tfse purposes of qwsliiative comparison. A 
< t a grade Photometries San%s CCD camera 21 equipped wilh a 

1300x1035 pixel Kodak KAM40O CCD ship was. added to the O.dx video port 
Ut s , oes of the sample 

x v * ^a of this embodiment of tho ^ 

Xv * ^e$> A^Bj^ermobrdrive system £2 was attached 

to the hxus knob of the microscope, This sbpper motor 22 was attached to the 
20 x — JH* Pamlum PC 23 also used to control the CCD camera 

ah enabhng tall automation of the acquisition of through-focus image 
sopuwxx - s ^ ^ s ^ ^ ^ 

1 recover pause images from the CCD images, thereby enahilng rati aafxmatlon 
>f IPC axxx ac ails t ; a a xi data pa 5 $ 

25 

In order la demonstrate that phase imaging using this Invent-on can accurately 
measure the phase structure of mcroseoplc samples it was accessary to have a 
sample with a wettxharaotehsed geemeay arid refractive index distribution. A 
3M FXSN-3224 optical fibre (a commercial fibre made by 3M) was used. 
30 xx c c measurements of the retractive index olsthaurwn N s \ oxinc 
atomic force necroscopy and commercial profiling teohnlaoes weru available 
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enabling aca ha - , •., , , , > , fi e kl 

Tnot - advantage of this bre was the 1 Tre v N 

effective ces ; an inner a maog ,a s m no w^eas^; 

fibres - mp ? # v the 

be v sr. I had to aocu^taiy .ir^age msltbns in 

refractive index rather than just two 

The optica liba >;so maged sfcle*on so as t, , v ^ , 

refractive index through all layers of the fibre structure. This was done by first 
stripping ?i:~ piastio sheath from the fibre by joskin 

using a commercial fibre atrpper to remove the plastic coating A small 
segment of fibre, typically a strgntf of approximately ens te two ceniimetres in 
>soopa a area - - < matching 

fluid and covered with a CWSmm thick cover glass, Anv hit , 
would introduce a xpunous tilt Into tlva recovered poase so two small sections of 
fibre, both of similar diameter to the sample, were placed parallel to and about 

am e air^s \ \ov s a\c i was as parallel fa „ 
-possible. 

images of the ®m were taKen oaing an Olympus 40x 0..7MA UplanApo 
objective, which meant that a 500x500 pixel imago conveniently spanned the 
whale width of ine l'-om : and - condenser was set at NA-0 2 The fibre's 
mbaohve index profiles were known fcr OTTSnm OAeTe laser) light so a 
825 ±1 Com band-pass Inteheranoo filter was Inserted into the Illumination 
sysfc * mt ap as profiles - a 

wavelength as close as passlbia to that for which data was available on the 
fibre. An Intensity image a! this sample in the plane of best focus and at ± Tom 
edher side of bast focus is shewn in Tome 8. alongside a phase image 

x Ki horn the « defoouse< ages using the phas< c wal a 
described above. Note that the fibre Is virtually Invisible la the livfocuc Image 
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^w w - \ o , - - ^ , w > n am to 

eg >ns of | are 

dearly visible in the phase image, 

i ; iau?6 9 shows a comparison of the measured and expected phase profiles with 
the uncertainties Indicated h the figure repressing eae standard deviation of 
the data along the length of the mm. This variation is thought to be enmahiv 
due to spatial variations in thickness of the cover g * a scope slkk 

As can he seer?, the recovered and expected phase profiles are In good 
agreement with one another, mm {he predicted 'profile lying within me error bars 
of the profile produced using the technique of this Invention. 

isisllu^ 

0 Pirates the appl rtiormAcg * < 
the three-dimensional imaging of objects throug! ov ,v v ( 
ton'iography techniques. 1 res is posrsole ss;ng the techniques ot this invention 
I ocou> > me eh we shsits ntmdueed by the object can ba oweth mewwod 

1 \ , atbns m the object, thus art m 
transform can be used to recover the three-dimer\ 

the projection data. Although theexperime\ ^ \. m, \ is in the 
optical regime, the same principles are equally applicable to X-ray, slow on and 
wutu s w 

For the purposes of collecting three dimensional data sets the same optical 
microscope described In the previous example is used with the addition of a 
rotation stage 24 for the purposes of rotating the sample within the con fines of 
the optical microscope Imaging area as shown in Figure 10. 'The rotation stage 
24 Is shown In greater detail in Figure 11. 
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The previously described arrangements inched a stepper meter ddve system 
22 attaehed to the parallel port of the same 133MM?. PmMum PC used to central 
the CCD camera 21 to drive the focus knob of the microscope A second 
steppe -note 25 was ) need 
5 system 24 for the purposes, of rotating the sample. This data acquisition system 
was linked to custom software written fa recover phase images from the CCD 
\pv thereby enabling full automation of trie mage wqulsaion and nets 
pmcemese sequence, Each data .set was collected using the same microscope 
as in Example 3 - an Qlympm BX-60 optica! microscope equipped with 3 set of 
io cover-slip corrected UpianApo objaollws and a universal condenser to provide 
KohPu iiiuudnatson, Jgitat Images vwwe cammed using a 12-bit Photometries 
SeuSys CCD samara equipped with a Kodak KAF-1400 1300x1035 pixel CCD 
chip on the 0.5.x video port of the mferoecepe, 

is To prepere the fibre sample 28 for Imaging the w - * * . 

a ^'oernern -J me end of a -section of fibre, typically .about one centimetre In 
\hmg the fibre m kopmpyl alcohol then using a commercial fibre 
atrspper to remove the plastic coating. The fibre was then cut into a small 
seamem of approximately one inch in length, with too ynstipped end then being 

20 slid into the end of a 26 gauge, 100mm sydhge needle 2? and fixed into position 
with a small amount of 8 minute A&Mtet. A mount 28 was osad to attach the 
needle 2? fe stepper motor 25. A pool of Index-matching fluid 29 surrounds the 
libre 20 as shown in Figure 11, with a microscope slide 30 affixed underneath 
the fibre using silicone grease arid a 0.15mm thick cover glass 31 placed evar 

25 the lap. 

Transmission wtenssly Images warn collected In the same way as described In 
Example 3 above using an Olympus 2 OA 0.4 hNA UMPian objective with the 
iden&s a MA-A ve images taken wen 
30 conve y spsnw , v of ihe fibre f 

precession. As the refractive ;ndey profile cf this fibre was Known for dba.dnm 



* laser) light, a 625 ± IQnm bandpass interference filter was inserted Into 
the illumination system to ensure that the recovered phase profiles were 
obtained o a wavelength as oboe no possible to that \ which data on the fibre 
was available. Baco pham fms 

either ssle of best foois, and data was coik ed boo , , - - 

through 180 elegies equally spaced in steps < area veen images A 
typical tomographic phase image fe shown in Figure 12. 

m ; ,--u-r data. rn. ti® fearo the mcork^trostod pmise images, ewm then 
processed into three-dimensional data seta using - xew> cocvsv-i-. 

m of the summation .of filtered back , „ w 

coda to perform me tomographic «K*5ii«t*uGtfcr> \ ft* 



to defermioe the looatlon of the notation axis and the data was digitally shifted so 
wied with the n«ddk co t sine 

^ • n iuctiort preset Fitting a serve to the phase profiles also 

nmeoo; - v v 4ata sets to be moved haci< t e. edde 
improved the quality of the ' ' ^ sled image. To s , w, 0 0 - - ^ - 
->lice through t\ 



Mote that all three regions of different refractive index are clearly resolved 
that these regions form conceetee cylinders, as is expected he tins sample 
m profile through the seams of this mroiedruotion Is shown In Figure 13 
sheet I on) alongside the Known refractive index distributor* for this fibre * solid 
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line), hue ea >n are very ok > those of 

the known profit which confirms the quantitative phase tomography technique. 




1. A mi - 

wave deid including the steps of 
(&) producing a r* r wsetvMm measure of the rata of change ef -a- c *\ 
5 of said radiation wave fold ov< n 

-As field, 

(b) produahg a repmsentaisve measure o* - , v vv 

wave field « elected surface; 

Co) transforming said measure of rate of eha ge < « , $ uce a 
to ? a transform repfeaen first 

i-tonrai transform representation a first Mar corresponding io the 
inversion of a first differential operator reflected in said me is m of 
rate of change of intensity to produce a first moo 

ts (d) applying an Inverse of said first integral transform to said first 
^odmsd itegral transfer?* epr* an 
* ^ - , presentation; 
(0) applying a correction 'based on aald measure ot inlet dfy ov« said 
selected surface to said unfransfor'med representation; 
a) - ; * erected untramforme< . 

second snfegrsi transform representation and applying to acid second 
aesenfafian a second filter corresponding to the 
E a second differential operator reflected in the corrected 



(91 



i produce a measure of 
selected piano. 
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transforms are produced using a Fourier transform, 
3, A method as claimed in claim 2 wherein said Feeder transform Is a 

Fast Former transform. 
5 4. A method as claimed in any one of claims 1 to 3 wherein said first 

and. seeofl 
opemtora, 

5. A method as claimed in any one of claims 1 to 4 > 

ssr ; 5 a some as mM second filter 
m a. A method as claimed in any one of claims 1 to e wherem said first 

filter includes se^ctsvely suppressing first h 



7. A method as claimed in any one of claims 1 to 5 wherein at least one 

- sa si and <■ end \ * \ das? correction fo 
n representative measure of Intensity, 

B.. A method as dinned In any one of claim* - t 

producing said mpresenfative measures of intensity and rate of 
cm Km ol merely over said selected surface fey producing 
representative measurements ^responding to intensify over at leas! 
20 two spaeed apart surfaces extending across the wave field . 

9.. A method as claimed in claim 8 wherein said selected surface is 

between two of said spaced apart surfaces. 
10. A method as claimed in claim 8 wherein said selected surface Is one 
mid 5 ed apart su ;e 
25 11. A method as claimed In any ana of claims S to 10 including the step 
of directly defecting representee measures of intensify over said 
spaced apart surfaces 
12, A method as claimed in any one of claims 8 to 10 including the step 
of producing said repn « y ove n east 

3Q - >a spaced apa t su fee es b> $ snnaec 
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13. A method as claimed in any one of claims 8 to 12 wherein said 
spaced apart surfaces are substantially parallel 

14. A m&hqd as claimed m dam 13 wherein said snacsd apart surface 
am s-. planar. 

s v c k mtho$ n < any «r*» of dale v » t< 14 m ? e * said 
reprosetru , ^ 

subtraction of mpteaantatfve tneasuremants of mtensiiy mapemmely 
madaaf legations over sakl spaced ar a ^ ^ 

16. A method as claimed m any one of claims 1 to 15 wherein said 
w representative measures of intensify and rate of change of intensity 

are obtained by sampling measurements at selected locations over 
said surface. 

17. A motel as claimed in claim 16 wnemn - s; . ; v 
messumn»rtls'a'r» made at locations defining a regular array over 

15. said surface, 

18s A method as claimed it* oiaim P , \*tu»** 

mm field propagates in a z-direotlon of artesian co-ordmat 
system: and further inducting tha step of producing an x componens 
and a y oomponant of phase separately. 

20 19. A method aa claimed in claim IS wherein aa>d first end sssi second 
filters have a component for product the x component at phase 
a a < i\ fiM producing the y component of phase of the 

form 

25 

o 2Lt&& 

wham k s .k, ere the Fourier variables conjugate to x and y: 
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a is a constant determined by oeise in the intensity measurement. 

A meiftod as ( iaimed n < ia m 19 N oknmx the so p ;! rj u \ ph ng s sld 
representative measure of rate of change of intensity by the negative 
cf average w« ve number < 
transformation. 
Amefotv ^ „ 

obtaining said representee measure of rate of change of intensity 
by obtaining a first ^res»MK measurement over a measurement 
surface across the wave field for relation of a first energy and 
draining a second teprmmdfoet measurement over said 
measurement surface for radiation of a second different energy. 
A method as claimed In claim 2 or claim 3 wherein at least one of said 
first filter and said secomi filter include a correction for aberrations in 

N v o x vasures of Intensity and rste of change of 
^ ^ v \ * ^ at least-one component (hpen4et\i on the 
aberration, coefficients of a system producing the repmoenmme 

A method as claimed in claim 22 wherein said radiation wave field 
propagates in a indirection of a Cartesian co-ordinate system and 
g the step of producing an x component a y 
component ct pease separately. 

A method as claimed In claim 23 wherein said first and said second 
filters have a component QJor producing the x component of phase 
arid a component Q y for prodaoiog the y component of phase doth of 
tli a term 

1 

v ^ k, an tlx Fourk sole > gafa ; < i 
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iVI} apply to said second integra! translouo representation a 
seoond •?•;:•-=• correspond inp to tho uwersioo of a second 
differential operate Sector a correct s c oed 
representation to produce a second niodsiod imepral transform 

(VII) apo's * " said 
seer x nvx » : ^ v ^ o produce 

measure of ptem of said radiation wave field across easd 
selected clone. 

ig 28. An apparatus as claimed fa oiafm 2? wherein said first and second 

, - > \ cims are produced using a rouxvn transform. 
29. An claimed in claim 23 wherein said Fourier transform is 

snsform. 

N - * v ^ * ^ - \ : \ x x x ^ said 
IS ! sand second Iffereoi re secoar rder differentia! 

oparatons, 

31 . An apparatus' as claimed in any on® of claims 2? to 30 wherem said 

frrsf filter is substantial^ th» same as said second filer, 
32 \n apparatus as claimed In any on© of claims 27 to 31 wmoein sasd 

20 x ^ v ^ s ^ " - \ \ V \ 

toe toot aitepno traostorni represents sore 
33. An apparatus as claimed in any one of claims 2? to 31 wherein at 

least one of said cod and second filters Includes a correction Ibr noise 
lo said representative measure of intension 
2:- 34 v N N - v 

ending to 

intensity over at least two spaced apart surfaces exfenceyp across the 
wave tield. 

35 kn apparatus as counted o 3 ndeol 

30 N - v s " as seaeod o. \> t ^ - v " ^ 



36, An x ^ ao claimed in claim 34 wherein said selected surface is 
one of said spaced opart surfaces. 

37, An apparatus as claimed in any one of claims 34 to 36 including 
detector moans positioned to directly detect rooooc-ndvOive measures 

5 ef Intensity over said spaced apod surfaces. 

S3 - * o < o N - % - s - i 

delec >r means to produce s 

over at least one of said spaced apart surfaces and imaging oceans to 
rosea ih a : r face onto the detector. 
w 33. An apparatus as claimed in any one of claims 34 to 33 orhe^n ssiu 
spaced apart surfaces a?e sebstanti * . * 
40, An apparatus as claimed in any one of claims 34 to 38 wherein said 

spaced apart surfaces are sue - % pi ^'ve 
41 An apparatus as claimed in any one of cke&m 34 1© ^ N * irs 
.33 means to produce i t n \\ m s c 

x d N ^ a st N s \i#&m measurements of Intensity 

respectively made at locations over said apaeed apart surfaces, 

reus as claimed in any one of alarms 27 to 4 1 wheresn said 
mean* to produce a representative measur 
so means to produce a representative measure of rata of chance of 

intensity sample at selected locations over said surface. 

43. ^n> ^ea ^ ,0xs\ ^ * ^ , 
at locations defining a regular array aver said surface. 

44. r - said 
25 v il^oo-; !rA v cam » . ■ Oartosian co« 

ordinate system and processing means produces an x component 

45. N , w >v> x x \ ' ' oo . ^ mm . m^. ^ s 
applies sasd first and said second filters nave a component O, for 

3tt producing the x component of phase and a component t for 

producing the y component ef phase of ihe farm 
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ohere lubto are the Fourier va da hies conjugate to v and y; 

a is 3 constant determent by noise in the intensity measurements. 

\n apparatus as claimed in claim 3? thereto as i representative 

nea^ure a? r&ie of change ef intensity - 

he average wave number o 4 



j ! wherein sale' 

represents* ^ . e of rate of cha x 
s a firat representative rnaaa^ < 
v - N *w field for radiation of , - 

obtaining a second representative measurement ever said 

An apparatus as claimed la eiaim 28 cr crs 29 ol worn at least one 
of said first niter and sale second filter include a correction for 
aberrations in sato representative measures of intensity and rate of 
change of imemtoy by including at ieast one component dependent 
an the aberration, coeltsclents of a system producing the 
represefbative measures. 

An apparatus as claimed m claim 43 wherein said radiation wave field 
propagates in a mtoreot;on of a Cartesian co-ordinate system and 
wherein an ,< concern- to and a y component of phase are 
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An apparatus as earned in claim 40 wherein said first and said 
second filters have a component hp for producing the x component of 
phase ana a component % for praduciop the y component of phase 
both of the form 



s- -r S ;: ) - 4T, ; Tp .d K Jpp;; 

! here k s dp are the Fourier variables conjugate to x and y; 
' is the average wave length of the radiation; 
Wi*t*-A*^} ® tne aberrated 'i$©fts8y 



61. \, 

x x set to a radiate % 

(b) nmesowg a repi^ntatlve measmo of s rate <*f change of 

wave field on the side ot the object remote from incident 



a measure m esensriy of rona .radiation 
v - acted surface; 
P.h r ov nsrg said measure - * to 

produce a first integral transform representation and applying to 
sa i - a first hue? 

:, espondmo ;. the inversion of a firs; Mar r 
rafieofed in sash measure of rate of change of intensity to 
groauce a wot won had « try - uwisform cm 
(e) applying an inverse of said first integral transform ta sees first 
^odsne v ^ \ s N N N N N s an 

untransformed repress nta lion , 
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of pciying erection based • eas sPy ove 

said selected - > i ^presentation 

(g) transforming the corseted urursnsbremed rennesorpaPon to 
produce a second integral transform representation and 
applying to said second in. , s a 

second filter corresponding to the inversion of a second 
differentia! operator reflected in the corrected ardransforrrmd 



(h) applying an inverse of said seco 

second modified integral tmnsform representation to produce 
measure of phase of 
selected plane* 

A metdoe ^ ^ v N x a Nv ^ v , v * v n s n 

o \ ui at ^ c*f hs gr ot ntens% 
over said selected surface by product epreserbaiive 

or spending to intensity o d leas 
apart surfaces extending acraaa the wave debt 
A mate 00 a sariss is 

between two or sani spaced apart sconces. 

A method as claimed In claim 52 wherein said selected surface Is one 
^ x v faces, 

A * any- one of da ms 52 s v a ^ 

snaced apart surfaces arc substantially parallel. 



representative Pleasure of rate of change of intensity is produced by 
subtraction of representative measurement? of intensity respectively 
mads at locations over sale spaaed apart surfaces. 
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58 A method as cla-nwd In claim SI including the step of producing said 

sakl 

selected , , 

correspe'\ K- t en&;ty ov< > ; ^tm 

s e I , xs oduc \ s coi d 

corresponding to intensity over said first strfe&e Cor tea changed 
distance he 

59. A method as claimed in claim 68 wherein said selected surfaces is 
10 - N m; . ! v * > rased apart surfaces, 

60. * method c c tmi i v y 

of direo% detecting reprosantaiiive measures of intensity over said 

81. A method as claimed in any one of claims 51 to 00 wherein said 
rs selected surface is spaced apart torn said tejeet m the ; roots w a 

propagation of said radiation, 
62 A method as Uaimed in any one of claims 51. to 81 wherein said 
a point source. 

63. A method as claimed In any one of claims §1 to 02 wherein said first 
20 and second integral transforms are produced using a Fourier 

ttan stores 

64. A method as claimed In any one of claims 51 to 63 wherein said 
Fourier transform is a Fast Fourier transform. 

65 do asperates fo laqln 

25 (a) a source to irradiate the object with a radiation wave field: 

(b) means to produce a representative measure of the rate of 
shaacs of Intensity of said radiation wave field over a selected 
s do* ^tending gear « fiek 

(c) means to produoo a representative measure of intensity of said 
30 radiation wave field over said selected su 

Cd) crocs as log means to sequentially 



(•) 5 - s - s Ot t V . 

produce a first integral transform representation; 
(ii) apod/ to said first Integral transform representation a first 

v\ 

operator reflected «n a ^ *\ of 

intensify to produce a first edified integral transform 



(M) t a to sslo first 

modified -Integral tram ^ , , N an 



* a - oijon bawd on so - ^ . , * 

said selected surface to said untransformed 
representation. 

(v) v * the corrected unpansforrned represantatioa to 

p; p oduot J se( | jsc ua? N ~ 

(vl) apply to said second integral transf< = epu - sntatii * a 
second filter coneap, ^ N 
differential operator reflected in the corrected 
ynfcansfem&d repmsenfation rc produce a second 
20 modified Integral transform representation; end 

p) apply an inverse of said second integral transform to earn 
second modified Integral transform representation to 
produce a measure of phase of sa;d radiation wove neks 
across said selected plane, 
as x s \ > s as 

mpresemtmve \ - ^ s corresponding to density ever at ieast 
two spaced apart surfaces extending asress the crave field. 
87, An apparatus as claimed is claim on wherein con selected surface is 

between two of said spaced apart surfaces. 
30 68. An apparatus as claimed in claim 06 wherein said selected surface is 
oneo s s iurfaces. 
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An " of claims 66 8 ( fini. 

detector means positioned to directly detect representative measures 
of intensity over said spaced apart surfaces. 

detector means to produce see rep; asem , moosme of intensity 
over at least one of said spaced apad sadaaas ana imaging means to 
Image thai surface onto the detector. 

An apparatus as claimed in an> a ^ m > x \ - , n\ 
spaced apart surfaces are sol 

An sippaisgfus m dallied fo -any one oi < s do i ; vdv > said 
representative measure of rate of ch,v ; — . \ ^ 

made at leaatioas over said spaced apart surfaces. 

An apparatus as claimed in .cte^ ^ m. s o, 

measures of ^tensity and rst< > Hie eosity cn 

* < b| producing firsi presenta uren 
corresponding to Intensity over a first surface extending across the 
wave field; means to change- fe distance between said source and 
sakl object, and means to produce second v ^ mau* 
masurc- ^ N s s \^ s , s\ s > ^ \x ^ 

between said o! 

An apparatus as claimed ;n claim 73 therein said aeiected sudace is 
one at said spaced apad surfaces. 

An apparatus as Osama in any one of claims da ta 72 including 
means to directly detecting representative measures of Intensity over 
said spaced apad sprfeaes. 

Ac apparatus as claimed in any one of claims 65 to 75 wherein said 
selected surface is spaaed apad from said object in ti 
; -op a; stlon of said radiation. 

An apparatus as claimed in any one of claims 65 to 76 wh- 
source Is substantially a paint source. 
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78. An apparatus as claimed ia any one of darns 66 to 7? wherein said 
first and second integral transforms ara prodaaad using a Fourier 

79. An apparatus as claims in claim 78 whmrem said Fourier transform is 
a Fast Foonet hailstorm 

80. A method ef phase amplifede imaging including x steps of 
(a) irradiating an object with a radiallon wave held; 

(p) foossslns -mm mm from tha object ihmegn an imaging system 
^ v v field 

a object; 

(c) producing a firs* representative mmmw of tnteratty 
distribution of radiation over s« at a first 

focus of the imaging system: 

(cJ) introducing a change in foeus of the image on said smaglne 
face tha 

(e) producing a second represemam/e of measure intensity 

distribution over so d imag face in 
(f; usmg said first ami second * - to 

produce a representative omasum of adenosy and a 
ve measure of rate of change of intensity in the 

dictation of radiation propagate; i surface 

extending across the wave field, 

m N w a tc 

produc* 1 - m • 

to ens - N filter 

corresponding to the inversion of a erst differential operator 
reflected in said meae.ee of rate of change of Intensity to 
a- aluoeaarsl m : era'? ep; meatalieo 

dm aopiymg an mvame of seal first integral Passman: to said ami 
modified miegmi transform representation to produce an 
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applying 8 correction based on said measure of eaenatty ever 
said selected surface to said unfraasformed representation; 
transforming the cerreatea aatransforniud representation to 
produce to said seeead Integral transform representation a 
second integral daastoae representation and applying a 
second filter ce 
differential operate 

r&^mmMScm to pw&tm a second modified integral transform 

applying an inverse et said second integral Oanstorrn to said 
second modified » \ . m 

measure of phase -of said ran. ^ 
) plane. 

81. A method as claimed in claim W wherein said radiaaon wave field 
has a numerieal ap&ium- smaller than the numerical aperture- of said 
lrr;agirn system. 

82 A method as daimed fct claim 80 or claim 81 said first focus of the 

n 19 ng system produces an intes image at the imaging surface 

and said second focus of the imaging system produces a shandy 

delecaseu image at the imaging surface. 
8.3. A method as claimed in any one of claims BO to 82 wherein sale 

irropyng surface is sudslaedtoy pisnar. 
04. A method as claimed In any oao cf claims 80 to 83 wherein free 

Imaging surface is an intensity detector, 
gc, A mcthed as claimed ;n any one of claims 80 to 84 wherein said 

imaging surface is said selected surface. 
88. A method as claimed in any one of claims 80 to 35 wherein said 

Integra! transform Is a Pearler transform. 
87, ^ . - s s x ^a v . x a s s ^ s ^ 

Fast Fourier transform: 



0) 
id) 
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88, s ^ ^ any em c N ) $3 id 

representative measure of rate of change of intensity is produced by 
subtraction of said first and second representative measurements of 
intensity* 

s • ; - 1 - , said 

representative measures of intveewy and rata of change of intensify 
are obfame ^ ^ ^ N ^ c ^ , 

said imaging surface. 

90. A method as claimed in claim 88 wherein said sampling 
so messommenie are made at locations v 

91. A method as claimed In claim 87 wherein see! radiation wave field 
propagates in a z-drrecSon of a cartesran co- radiate system and 
further including the step of producing an x component and a y 

15 component of phase separately, 

92. A method as claimed in claim 01 whe - 

fitex nave a component 0^ for producing rhe ,< component of passe 
and ■> component Q } , for producing fee y component of phase of the 
term 



where leeA are the Fourier variables conjugate to x and y; 
25 eea ie sun me s w 

93. ^ v e> x < ding the she 

representative measure of rats of change of Intensity by the negative 



of the average wave number of the radiation besore said integral 
traraaormaaon, 

S4, An apparatus for phase amplitude imaging of en apieci including 

a radiation wave held seume !o eradiate sa:d ooieat: 
an imaging system to locus radiation Porn said object to an 

imaging surface < x 

object; 

means to produce n eprese P >t vc m< a« ae of Kii >tk n 
^ \- * ovos - mirfaea; 

said imaging system including so a - eneraele means to 
adjust said fecus of said radiation to sad waning surface to at least 
at a first focus and a second focas; 

ens to, 

(s) produce a representative measure of nt^v; x ami a 
representative me&urs of rate of no ^ , ^ a ^ a - t - 
dinsetfcn. of radiation propagation over a selected surface 
extending aeross the wave field front representative 
measures of radial ntonsify ova* as a mage surface at 
- v^s and- said seoemd focus; 
go sanstorm said measure of rate of change of intensity to 

proouoe a nrst Integra; transform representation; 
gfu apply ta seat lea?, iemgral transanal representation a fast fiitar 
corresponding to the inversion of a first dibomntiai operator 
reflected In said measure of rate of change of intensay to 
^ ^ . s , ^entotson; 
(iv) apply an inverse of said first integral transform to said first 

p- ode? e an 

untfansfonaod representation; 
(v> apply a correction based on said measure of intensity over 
said selected surface to va:0 Cornea representation; 
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(yi) transform this corrected unimnsforrnod representation to 

odor >r agral tram on m^tadcn 

(Al) apply to said second integral transform representation a 
second filter corresponding to the inversion of a second 
. \l operator fefteed In the cor reeled mdranrdormed 
e N ^ it eon ic i kx 
transform represent; . s ~ one 

apply as m so of $a*f secv erne 
second mod&fed. Integral Iran )mi men la; > prota* 
a measure of phase of said mm * . - ^ - 

selected plane, 

95. An apparatus as claimed In claim 94 wherein said mention wave field 
has am ; ml aperture smaller than d\ i < 
m system, 

00 An apparatus as claimed m elate 94 or claim 95 sard era? feces of the 

Imaging syciom produces m infOCUS Image al do eacino cerfar.:e 
and said second pern of the imaging system produces a sightly 
defcensed image at die imaging surface. 

ens as claimed In any one of claims 94 to 96 wherein said 
anamne sunace a; sudsraneaAy planar. 

98. An apparatus as claimed m any one of claims 9a la a? wherein me 
imaging acrface is on intensity detector.. 

no An apparatus as claimed in any one of claims 94 to 98 whenon said 

imaging surface is said selected surface. 

100. An apparatus as claimed in any one of claims 94 to 99 wherein said 
Integral transform Is a Fourier transform. 

101. N N s 
Is a Fast Fourier transform. 

102. An apparatus as claimed in any ^ \ of claims 94 to 0 \ whonm> sasd 

On^I\V ^O ^ \ N V V v S ^ X . \ ^ 



subtraction of said test and second representative measurements of 
intensity, 

103- J V - - < } \S s 

*epreseev 

are obtained by sampling measurements ai. seiested location ever 
sssd Imagine scenes: 
104 An apparatus as claimed in elate 103 .wherein sed sampling 
measurements am made at locations defining a regular army over 
east imaging surface. 

105, An apparatus as claimed ift claim 10? * 

field propagates in a z -•direction of a cartesian co-ordinate system and 
c t is Step of producing an > component and a y 
*n ^ \ , s , w> separately 

106, An apparatus as claimed In daim 105 whareln said first and said 

, m a component njor pr< 
*,la eooipene a 
phase of the form 



X ~ u i 



where k,og are the Fourier variables conjugate to x and y; 
i.? ;u s cunu-evv. oei.es sees st ey sx^e :vsry ^ssas^^vs^e;;^, 

2$ N \ method as claim s — a of multiplying. 

said representative measure o? rate of change of mteusuy by the 

\ ve o € average wave w c \ 
; r see ra I trs n storm anon . 



A computer program for prooessing a representative measure of the 
rate of change of intensity of a radiation wave field over a selected 
surface extending o, across the wave field arte a 

t^rmmmm measure of intensity of said naslaimn wave held over 
said seeded surface, said program including 

:) « : i a ts! > sing said measi re of one oi diange of intensity to 
produce a ♦ , 
first integ * 

wersion of a first mt * - sa 
rate of change of mtmmty to produce a first modified Integral 
trans tort";"; ropro s en ia s on; 

«od» tor m¥m m inverse ©f said first Integra! transform fo said first 
modified integral transform egret arson : > \, an 
y ntrao $ % sedation j 

code for applying a eorr^tlort based on said measure of intensity 
over sale seleotec surface to said nntransformee representation:; 
code form ^ dx* corrected untianisformed re? osentatlon to 
prod ace a second integral transform representation and applying to 
ei.no ;n c .mi integral Iransfdrm repreaerstatlon a second filter 
corresponding to the inversion of a second dlffe t 
reflected corrected untransformeo r ^nn > o piedneo : 
second meditate integral fransiorm represerdahon; 
code for agpiying an inverse of said second Integra! transform io said 
second modified integral transform representation. 
A computer program stored on computer readahle storage media for 

ste o diange of 

intensity of a radiation wave item! over a selected surface extending 
generally across the wave field and a representative measure of 
intensity of sasd radiator; wave field ever said selected surface, said 
program including: 



sods for transforming said measure of mte of change of intensity to 
produce 3 first integral transform mpressntador; and applying | 0 sard 

hid atsform representation < jo 
Inversion of a first differentia? operator reflected in said measure of 
rate s >ge of N« >S tj b p < 
f ran storm re a re ss n ta i i on : 

59 an mmm of m6 first lidegmi transform m said first 
•noatfea itegra , $fom BpreK stson c produce an 

q a cormctkm based on vk >m \ , si v ^ 

code for tmsiaforming the cone 



corresponding to the inversion of a second differentia! operator 

second modified integral transform repmsantatfem; 
cade far applying an inverse of said second integral transform to said 
^ ^ ^ oor«enialte. 
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(c) negative etefee 
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?'fe) apertere pkr.se phase 




(1) recovered phase 
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